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A Toy Model of Sea Ice growth 


A.S. Thorndike 

My purpose here is to present a simplified treatment of the 
growth of sea ice. By ignoring many details, it is possible to obtain 
several results that help to clarify the ways in which the sea ice 
raver win respond to climate change. Three models are discussed. 
TTie first deals with the growth of sea ice during the cold season. 

^ ZZl/ sc 7 7 ? e 7 cle of growth 3X1(1 meItIng for Penial 

i^Vh ‘ i u CndS thC SeC ° nd tD aCCOUnt for the Possibil- 
ify that the ice melts away entirely in the summer. In each case, the 

o jective is to understand what physical processes are most impor- 
tant, what ice properties determine the ice behavior, and to which 
climate variables the system is most sensitive. 

Climate 

Y y T2 nt ° 3 C ° ld 311(1 3 W3rm season - each of duration 

tion K f g L C u Sea ?° n ’ the 3Ver3ge downwe Hing longwave radia- 
on is f, wc which equals 180 W/ m 2 in today’s climate. In the warm 

season, the downwelling radiation is f lww = 270 W/ m 2, and there is 
shortwave radiation f sw = 200 W/ m 2. fThese values are from Maykut 
and Untersteiner, 1971.) A fourth climate variable is the heat sup- 
plied to the ,ce from the ocean, f w . Its value is not knowTLZ 

th U ronih V n ar t ih ^ SPaCC and timC - WC Wl11 regard U as instant 
n wf 2 U ^ yC3r ’ 311(1 expect if to have a val ^e in the range 0 to 
/ m m the central Arctic. The four parameters f, w „ f, f 
and f w specify the climate. Iww ’ sw ’ 

To examine the response of the ice to changes in climate, we will 
ve n mind perturbations to the longwave fluxes, such as might 
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accompany tropospheric warming or changes in atmospheric 
composition: 

flwc ^lwc + ® 
f]ww flww + ° 


Sea Ice Parameters 

The ice will be described by two variables, its thickness h and 
surface temperature T. Other ice properties and their values used 

this paper Eire: 


longwave emissivity 
thermal conductivity 
specific heat 
latent heat 
albedo 


e 

k 


a 


2 W/m/K 

2 x 10 6 J/m 3 /K 

3 x 10 8 J/m 3 
0.65. 


In taking these parameters to be constant, we ignore the strong 
dependence some have on temperature and salinity in effect we 
renlace the thermal properties of sea ice with those of fresh ice. In 
so doing, we give up the possibility of resolving the sea ice beha o 
within a few degrees of the freezing point. 

Ice Growth 

The black body radiation from the ice surface can be expresse as 

iz: 

a is the Stefan-Boltzmann constant, 5.7 x 10 8 J/m /s/K , 

27 1 .2 K is the freezing point of sea water. 

In the cold season, and in the absence of upward heat flux from 
beL J the radiation from the ice surface must balance the down- 
welling longwave radiation. This implies 

T = -(A - fiwc)/B = -D/B = -29°C (1) 

, i^a convenient scale for the radiative 

toesTS the net radlLuon balance over a surface at the freezing 
point. Similarly, the quantity -D/B is a natural ^ ^e surface 
temperature. It is the minimum temperature the ice surface can 

31 Fot Ice of thickness h, having a linear temperature profile, heat is 
conducted upward at the rate -kT/h. and this flux must be Includ 
in the surface energy balance: 
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A + BT = f lwc - kT/h (2) 

It follows that the surface temperature and the Ice thickness are 
related as 

T = - Dh/(k + Bh), and h = - kT/(D + BT) (3) 

These relations hold when the ice is growing, because then the tem- 
perature profile through the ice is approximately linear, and the 
expression for the conductive heat flux is justified. 

At the bottom of the ice, the rate at which heat can be conducted 
through the ice determines the rate of ice growth: 

L dh/dt = - kT/h - f w ( 4 ) 

After substituting for T, the equation can be integrated to give the 
relationship between the thickness h, after time t, and the initial 
thickness h Q : 


tB/L + kDf- 2 log 


k(D-f w )-f w Bh 

k(D-f w )-f w Bho 


+ B(h-ho)f-i =0 


(5) 


In the special case ho - 0, the logarithm can be expanded in the 
form log (1 — P) =r — p — 02/2, provided f w « D, leading to 


L dh/dt = - kT/h -f w (6) 

This result establishes k/B as a natural scale of ice thickness. B 
determines how much the surface must cool to maintain radiative 
equilibrium. Together with k, this fixes the conductive heat flux, and 
therefore the ice growth. The value of k/B is about 0.4 m. It also 
appears from this result that kL/BD is a natural time scale, having 
the value of about 1 1 days. It is worth noting that the ocean heat 
flux decreases the thickness through the linear factor (1 - f w /D). 
With D a 140 W/m 2 and f w a 1 W/m 2 the ocean heat flux has a 
small effect on the ice produced in one year, but we will see below 
that it can have a large effect on the equilibrium thickness. 

It can also be noted that a steady state solution to Equations (2) 
and (4) is possible. The ice must grow to a thickness where the con- 
ductive heat flux equals f w . This occurs when h = (k/B)(D/f w - 1), 
which is about 50 m. However, the time required to approach this 
equilibrium is much longer than a single growing season, so the 
model needs to be extended to account for the seasonal cycle of 
growing and melting (see “A Model of the Perennial Ice Cover"). 

For very thin ice, t <sc kL/BD =11 days, we have the further 
approximation that 


h = (D-f w )t/L 


(small t) 


(7) 
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In this case, the conduction through the ice is not important. The 
ice growth is determined completely by the net heat balance for the 
slab, and the latent heat, and therefore grows at a constant rate. 

For large t, t » kL/BD, we obtain 

h = (k/B)(l - f w /D) V"2BDt/kL (large t) (8) 

which shows that h grows as t 1 / 2 . The thicker the ice gets, the 
slower it grows. Here the growth is regulated by the conduction of 
heat through the slab. With time, the ice both cools and thickens, 
but the combined effect is that the temperature gradient decreases, 
which controls the rate at which heat can be removed from the bot- 
tom surface, and thus the rate of ice growth. 

The ice thickness depends on the climate through the variables D 
= A - f lwc and f w . It is interesting to note that in the coldest climate 
possible, f lwc and f w = 0, the ice would only grow to h = 
(k/B)(2kAY/BL) 1/ ' 2 = 4 m in a single growing season, where Y - one- 

half year. , 

Now consider the case ho > 0, as for ice which has survived the 

summer. During the subsequent growing season, the ice will reach a 
thickness h after time t. In the case where the new growth h - ho is 
small compared to ho, we obtain 


h-hc 


kD 


k + Bho 


-f I— 

1 W 1 - 


(h-h 0 « ho) 


(9) 


This implies that the ice cannot grow thicker than (k/B)(D/f w - 1), a 
quantity that is sensitive to the ocean heat flux f w . Compare this 
result with Equation (25), below, which accounts for the annual 
cycle of ice growth and melting. 

Equation (5) can be solved for h numerically. Thus we can 
express 

h = H(hot) ( 10 ) 

though we cannot write the function H explicitly. 

The assumption of a linear temperature profile through the ice 
has been essential in obtaining these results. As the ice grows, it 
also cools, according to Equation (3). There is an internal inconsis- 
tency here. An element of ice cannot cool if the temperature profile 
is linear there. The only way to avoid this kind of inconsistency is to 
formulate the ice growth as a heat diffusion problem, as Maykut and 
Untersteiner did. To estimate how large an error the assumption 
introduces, note that the model properly accounts for the heat 
released by the new ice growth, Ldh, but it does not account for the 
heat which must be removed to cool the Ice, dQ, where Q = chT/2. 
For ice that grows one meter in the winter, Ldh is 3 x 10 8 J/m 2 . The 
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heat storage term is (2 x 10 6 J/m/K)(l m)(-21 K)/2 = 2 x 10? j/ m 2. 
So we may expect errors of about 5 to 10% to arise because we have 
not properly accounted for the heat storage (although in the next 
section it is shown that a large part of the heat storage can easily be 
accounted for). 

A Model of the Perennial Ice Cover 

Since the ice is characterized by its thickness and surface tempera- 
ture, it is useful to plot its behavior as a trajectory in (h,T) space. 
Figure 1 displays the results of the Maykut and Untersteiner model 
this way. We will approximate the annual cycle as follows. Imagine a 
slab of ice at the end of the melt season, having thickness h and 
temperature T = 0 throughout. When the weather turns cold, the ice 
cools until its surface temperature reaches the temperature appro- 
priate for its thickness. Equation (3). Then it grows for the remain- 



Figure 1 . The equilibrium cycle for the standard case of the Maykut and 
Untersteiner model, plotted in (h,T) space. 
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ing part of the cold season, according to the ice growth model of the 
previous section. In the warm season, the ice first warms up to the 
freezing point, and then melts. The cooling and warming processes 
take place with no change in thickness. They account for the stor- 
age of heat by the ice. The equations governing this four-process 
cycle of cooling, growing, warming, and melting are 


Cooling 

D = A — f[ wc 
T = -Dh/(k + Bh) 

Q = -chT/2 

x c = Q/(A + BT/2 — fj wc — fw ) 


(ID 


Growing 

Tg=Y-T c 

h g =H(h,T g ) 

T g=' Dh g/ ( k + Bh g) 


( 12 ) 

(Equations 5 and 10) 


Warming 

Q = -chgT g /2 

D = fjww + (1 - ot)f sw + fw - A - BTg/2 
t w = Q/D 

Melting 

T m = Y-X w 

h = hg -x w [f] ww + (1 - cx)f S w + ” A]/L 


(13) 


(14) 


where x is duration with subscripts c, g, w, and m indicating the 
four processes. 

A short program can be written to do these calculations. The 
results (see Table 1) show that for the present climate, the ice is 
attracted to a periodic orbit, as sketched in Figure 2. The ice is 
saved from the two alternatives, melting completely or growing with- 
out bound, by the negative feedback provided by the ice thickness. If 
the ice is too thin in the fall, it will grow more in the winter than it 
melts in the summer. If it is too thick in the fall, it grows less than it 
melts. In either case, it approaches the equilibrium orbit. 

For other climates, the orbit changes. For large enough perturba- 
tions to the climate, the equilibrium cycle cannot be maintained. If 
the climate is too warm, the ice melts completely in the summer; the 
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Table 1: Results of model simulations 


Perennial Ice Cover 

8 

(W/m 2 ) 

fw 

(W/m 2 ) 

(d) 


^ g 
(m) 


k) 

'w 

(d) 

(mj j 

-15 

1 

43 


9.4 


-32 

44 

9 1 

-10 

1 

21 


5.3 


-30 

23 

4.8 

0 

I 

7 


2.9 


-26 

1 1 

1.9 

5 

1 

3 


2.5 


-25 

8 

1 i 

10 

1 

1 


2.2 


-24 

7 

0.6 

1 5 

1 

0.1 

2.0 


-22 

6 

0.2 

0 

2 

5 


2.7 


-26 

10 

1.6 

0 

5 

3 


2.3 


-26 

8 

1 0 

0 

10 

0.5 


1.8 


-25 

6 

0.3 

Seasonal Ice Cover 








8 

(W/m 2 ) 

fw 

(W/m 2 ) 

^ml 

m 

(d) 


hg 

(m) 

& 

(d) 

(d) 

20 

1 

2 

37 


1.7 

-21 

5 

36 

30 

1 

6 

117 


1.0 

-16 

2 

1 15 

37 

1 

9 

170 


0.3 

-9 

0.4 

165 

20 

5 

4 

82 


1.3 

-19 

3 

82 

20 

10 

6 

119 


0.9 

-18 

2 

119 

h „ and T„ are the ice thickness and temperature at the end of the growing season. h m is the thickness 
at the end of the melting season. t c , V t ml . and t h are the durations of the processes which cool and 

1CC ' C °° r a / ld h f at the mixed la y er - 8 ls the perturbation to the longwave radiation fluxes. 

f w is the upward flux of heat from the ocean. T ml is the temperature of model's upper ocean layer. 


model is patched up to treat this case below. If the climate is too 
cold, the ice grows without bound. The behavior of the model as a 
function of the perturbation 5 is sketched in Figure 3. The longwave 
fluxes can vary up to about ± 20 W/m 2 from the present climate 
and still support an annual equilibrium cycle. 


A Seasonal Ice Model 

If the ice melts away before the end of the summer, the positive 
heat balance at the surface must cause the upper ocean to warm. 
The model of the preceding section can be modified to account for 
this heat by including an upper ocean layer. The thickness of the 
layer h m | is assumed to be fixed at 50 m, and its temperature is 
allowed to change. The layer is assumed to be well mixed, so that a 
single temperature describes its state. Begin the annual cycle at the 
end of summer with the mixed layer having temperature T mI . We 
calculate how long it takes for the mixed layer to cool, and then how 
much ice will grow during the remainder of the cold season. During 
the warm season, the ice must warm up to the freezing point, and 
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Figure 2. The equilibrium cycle for the perennial ice model, under the 
present climate; C is cooling, G is growing, W is warming, M is melting. 


then melt, requiring times x w and x m . The remainder of the warm 
season is spent heating the mixed layer. The model equations are. 

Cooling the mixed layer: 

Q = c w h ml T ml (15) 

x m ] = Q/( A + BT m j /2 - fi wc - f w ) 


Ice growth: 

Tg = Y — t m | 

h = H(0.t g ) 

T = -(A-f lwc )h/(k + Bh) 


( 16 ) 
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Figure 3. The maximum and minimum thicknesses in the annual cycle for 
the perennial and seasonal ice models, as Junctions of the perturbation S to 
the longwauejluxes. 


Warming: 

Q = -chT/2 (i 7 ) 

D = f iww + (1 - a)f sw + f w - A - BT/2 

t w =Q/D 

Melting: 

~ hL/[fiww + (1 ~ °Ofsw + fw — A] (18) 

Heating the mixed layer: 

T h — Y ~ ~ T m 

^ = flww + (l “ a w )f sw + f w — A 
T ml = (°/ B )[l - exp(- Bx h /c w h ml )] 


( 19 ) 
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For sea water the specific heat is c w - 4 x 10 6 J/m /K and the 
albedo is taken to be 0.2. The last equation is the result of integral- 
ing the heat balance for the mixed layer. 

d(c w h ml T ml )/dt = D-BT m i 

For positive perturbations to the longwave fluxes in the range 17 
W/m 2 < 5 < 40 W/m 2 , a seasonal ice pack is possible (see Figure 3). 
For example, with 5 = 30 W/m 2 , the mixed layer reaches 6°C. It 
requires 1 17 days to cool off in the fall. In the remaining 65 days oi 
cold weather, the ice grows to 1 m, reaching a surface temperature 
of-16°C. In the warm season, it takes 2 days to warm the ice to t e 
freezing point, 65 days to melt it, and 115 days to warm the mixed 

^For climates warmer than 5 = 40 W/m 2 , the mixed layer gets so 
warm that it cannot cool completely in half a year, so ice never forms. 


The Annual Energy Balance 

The coldest surface temperature and the thickest ice occur at the 
end of the freezing season. A simple energy argument leads to esti- 
mates of these quantities. Since the temperature varies between T - 
0 and T = T fi during the cold season and during the warm season, 
the radiated energy for the entire year must be approximately 

2Y(a + BTg/2) 

In an equilibrium cycle, this must equal the energy reaching the ice, 
which is 

Y (fj wc + fw ) + Y [flww + (1 “ a ) f sw + ] 

Equating these, and neglecting f w for the moment, leads to 

T g =(D m -D g )/B 


( 21 ) 


in which D g and D m are the net radiation balances over an ice sur- 
ace at the freezing point during the growing and melting seasons. 
Jsing this expression in the surface heat balance. Equation 
fives the maximum thickness as 

( 22 ) 


h g = (k/B)(D g -D m )/D m 


With D m = W + (1 - a)f sw - A = 20 W/m 2 , and D g = A - fj wc = 140 
W/m 2 , we obtain T g = -26°C, and h g = 2.6 m. The expression for h g 
underlines the sensitivity to the heat balance in the melting season. 
As D m approaches zero, the ice grows without bound. 

The simulations in “A Model of the Perennial Ice Cover” showed 
that the perennial ice cover gives way to a seasonal ice cover or a 
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rr per j urbation of 18 W/m 2 to the longwave fluxes. To see why 

^eVnnl ?/ a , COnditIOn on °g D m which just causes the 

ce to melt completely at the end of the melting season. If h = 0 at 

e eginning of the cold season, the thickness at the end of the 
season will be, Equation (6), 


h = (k/B) (l + 2BYD g /kL) V2 • 


(23) 

In an equilibrium cycle, this must equal the ice melted in the warm 

season which is approximately YD m /L. Combining these results 
gives the constraint 

(! + BYD m /L) 2 = l + 2BYD g /kL (24) 

Pl ° tte n in f igUrC 4 ' F ° r Climates havln g ( D g. D m ) below this 
n - Jo 1CC S l melt aWay com P le tely during the summer. If we set 
m +6 and D g - 140 + 5 the equation requires 5=17 W/m 2 in 
good agreement with the simulations. 


Sensitivity to the Ocean Heat Flux 

One result of the Maykut and Untersteiner simulations is a strong 
dependence on the ocean heat flux. They showed that Increasing f 
from 0 to 7 W/ m 2 caused the equilibrium thickness to decrease 
from about 6 m to zero. When one considers how poorly the ocean 
eat flux is known, and the uncertainties in the other heat fluxes 
this sensitivity is unsettling. 

The present models have been run using different values of f 
with results that support the Maykut and Untersteiner results. In 
particular, holding other climate variables fixed, f w = 1 W/m 2 nro- 

w/^ t °\ . t ^ ckness 1-9 m at the end of summer. With f w = 12 

W/m 2 that thickness is reduced to 0. 1 m. 

We can appreciate how this works by restoring f to the expres- 
sions for Tg and h g . Equations (21) and (22): 


T g (^m D g+2f w )/B (25) 

h = k (D g -D m +2f w ) 
g B (D m +2f w ) 


From the first expression, we see that the surface temperature 
responds to the ocean heat flux exactly as it responds to the radia- 
tive fluxes. If an increment of heat A is added to the ice, in any way 
the ice must respond by warming its surface A/B to maintain the 
overall energy balance. However, the response of the ice thickness is 
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d) m : Mgt ‘Radiation during Melt (‘W/m 2 ) 


Ftnure 4. The condition on the net radiation balances in the cooli ^ g a ^ d 
mMtna seasons D a and D m . which causes the ice just to vanish at the end 

hne represents climates cn which the longwave 

fluxes differ by 5 from the present climate. 


nore subtle. An Increase to f w decreases the numerator ^and 
ncreases the denominator, and therefore has a larger e ec on 
hickness than one might expect. 


The Time Scale for the Ice Response 

By combining the ice grown in the cold season and the melted 
in the warm, we obtain the ice thickness after one year, beginn g 
with initial thickness h n . 


hn + l= F ( h n) = h n +YL_1 


kD 


k + Bh n 


2f w -D r 


(26) 


Now the time scale for the Ice to reach Its equilibrium cycle Is 
r = h n -h eq /(h n -h„„) (years) 


( 27 ) 


ft.s. ‘Tfwrtufife 


237 


If we also write 

hn+i — h e q +(h n — h e qjdF/dt 
and substitute for h n+1 in (27), we find 

r = (l-dF/dh n ) _1 

where the derivative is evaluated at the equilibrium thickness 


h eq =(k/B)| 

This gives the time scale 


D 


V +2f w 


-1 


T = (kL/YB) 


D 


(D m+ 2f w ) : 


= 3 yr 


(28) 

(29) 

(30) 

(31) 


Discussion 

The essential mechanism at work is that the ice adjusts its sur- 
face temperature to maintain the energy balance at the surface. In 
the present climate, the annual energy balance over a surface at the 
melting point is D m - D g = -120 W/m 2 , meaning that a surface at T 
= 0°C would lose 120 W/m 2 more than it receives. The surface 
adjusts by cooling to about -26°C. 

The second important idea is that during ice growth, the ice 
thickness is related to its surface temperature; the colder the tem- 
perature, the thicker the ice. Therefore we can assess the response 
of the ice thickness to a climate perturbation, in the form of a 
change dE in the energy reaching the surface, by evaluating 

dT/dE and 

dh/dE = dh/dT • dT/dE = -[(k + Bh) 2 /DBk] (32) 

For h = 3 m, the sensitivity is about -0.2 m for 1 W/m 2 increase in 
incoming energy. Note that the sensitivity is a strong function of 
thickness. 

Thus the thickness is rather sensitive to the energy fluxes. Pertur- 
bations of 10 W/m 2 can mean several meters of ice thickness. Never- 
theless, the system can adjust to perturbations of this magnitude. 

However, for larger fluctuations, the system cannot adjust. A 
cooler climate in which the longwave fluxes were 20 W /m 2 less than 
present (about 5°C of tropospheric cooling) would allow the ice to 
grow without bound. A positive perturbation of the same magnitude 
would cause the ice to melt away completely in the summer, so that 
the entire Arctic would support only a seasonal ice cover. With a 
positive perturbation of about 40 W/m 2 , enough heat Is stored in 
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the upper ocean during the warm season to prevent ice from form- 
ing during the winter. 

Of course, with a model that ignores many physical processes 
and treats others badly, one cannot have great confidence in the 
numerical results. Still, I expect the relationships between quanti- 
ties in these models to have a qualitative validity, in the sense that 
in a more careful model, the signs and powers in the expressions 
would survive, while different numerical factors, of order unity, 
might appear. 

These ideas have implications for climate models that attempt to 
simulate the Arctic ice cover. First, the level of detail in the peren- 
nial and seasonal models discussed here may be about right for use 
in a climate model. The models surely allow an interactive ice cover 
without a heavy computational burden. A second point is that we 
have found limits on errors that can be tolerated in the energy 
fluxes supplied by the climate model to the ice. Fluxes in error by ± 
10 W/m 2 will not destroy the ice, but errors of ± 20 W/m 2 will lead 
to a qualitatively wrong ice pack. Finally, given the sensitivity of the 
ice thickness and temperature to perturbations on the order of 10 
W /m 2 , and the likelihood that climate models will have this much 
uncertainty for some time to come, it is not realistic to expect a cli- 
mate model to reproduce the observed ice thickness, or to predict 
how much the thickness will change. It is enough to ask that a cli- 
mate model distinguish between a perennial ice cover and a sea- 
sonal ice cover, or between an ice cover and no ice at all. 

Of course, the actual Arctic ice cover is too complicated for a toy 
model. There is in fact a great range of thicknesses, all growing and 
melting at different rates, and all responding to mechanical as well 
as thermal forcing. I do not want to leave the impression that these 
processes are unimportant. 
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